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Abstract. We derive a spin-orbital Hamiltonian for a triangular lattice of Cg orbital 
degenerate (Ni^+) transition metal ions interacting via 90° superexchange involving 
(0^~) anions, taking into account the on-site Coulomb interactions on both the anions 
and the transition metal ions. The derived interactions in the spin-orbital model are 
strongly frustrated, with the strongest orbital interactions selecting different orbitals 
for pairs of Ni ions along the three different lattice directions. In the orbital ordered 
phase, favoured in mean field theory, the spin-orbital interaction can play an important 
role by breaking the C/(l) symmetry generated by the much stronger orbital interaction 
and restoring the threefold symmetry of the lattice. As a result the effective magnetic 
exchange is non-uniform and includes both ferromagnetic and antiferromagnetic spin 
interactions. Since ferromagnetic interactions still dominate, this offers yet insufficient 
explanation for the absence of magnetic order and the low-temperature behaviour of 
the magnetic susceptibility of stoichiometric LiNi02. The scenario proposed to explain 
the observed difference in the physical properties of LiNi02 and NaNi02 includes 
small covalency of Ni-O-Li-O-Ni bonds inducing weaker interplane superexchange in 
LiNi02, insufficient to stabilize orbital long-range order in the presence of stronger 
intraplane competition between superexchange and Jahn- Teller coupling. 



1. Introduction 



The low-temperature magnetic behaviour of LiNi02 has remained puzzling ever since its 
peculiar properties were discovered p. For no apparent reason, LiNi02 does not show 
magnetic order nor a cooperative Jahn- Teller effect down to the lowest temperatures, 
which is very different from the conventional behaviour observed in its sister compound 
NaNiOa. 
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Figure 1. Structure of LiNi02: (a) Fragment of 3D crystal structure, with Ni ions 
shown by filled circles, Li ions by gray circles, and O ions by open circles, (b) Nickel 
layer with adjacent oxygen layers. Each Ni^+ ion (black circles) is surrounded by six 
oxygens (open circles). The directions of the Ni— Ni bonds are labelled as XY , YZ 
and ZX, corresponding in each case to the plane spanned by the connecting Ni— O 
bonds. 



Structurally the two systems represent an interesting special category within the 
class of correlated transition metal (TM) oxides. LiNi02 has a layered structure [see 
figure m^a)]: it is rhombohedral, consisting of successive (111) planes occupied by Li+, 
0^~, Ni^"*", and ions. Thus the Ni^"*" ions are on a triangular lattice, with each 
direct Ni-Ni bond lying along the diagonal of a nearly square Ni-O-Ni-0 plaquette, 
the Ni-O-Ni bonds being close to 90 degrees. This is distinct from the more commonly 
encountered situation where the bond between two transition metal ions through the 
ligand ion connecting them is close to linear (180 degrees), as e.g. in the perovskites. 
As pointed out by Mostovoy and Khomskii this difference should have important 
consequences for the orbital and magnetic superexchange (SE) in LiNi02, since the 
SE within a Ni plane with Wi^+{tlgel) ions would originate predominantly from virtual 
charge transfer excitations e^2p^e^ ^ e^2p^e^ ^ e^2p^e^ along the 90 degrees Ni-O-Ni 
bonds. 

Over the years, a number of experiments (magnetic susceptibility, ESR, NMR, 
neutron scattering) have been performed on LiNi02 samples of varying stoichiometry 
j3]-JT]. From these data one has concluded that the presence of excess Ni ions in the 
lithium layers introduces extra ferromagnetic (FM) coupling between the nickel layers. 
In addition, for the samples closest to perfect stoichiometry a positive Curie- Weiss 
temperature was found, indicating that the in-plane exchange is also FM. This is not in 
accordance with the description given initially by Hirakawa a/ jHI , namely that LiNi02 
would be a triangular lattice antiferromagnet (TALAF). Actually, this assumption was 
the original motivation for performing magnetic measurements on LiNi02, since the 
TALAF for spin = 1/2 is a frustrated system and the ground state might be some 
kind of quantum liquid |12] instead of the classical 120° rotated spin arrangement. 

At first sight the FM nature of the in-plane correlations is not surprising when one 
looks at the three-dimensional (3D) crystal structure (figure Q). As the neighbouring 
Ni^"*" ions are connected via two Ni-O-Ni bridges, one has the case of 90° SE, which the 



Orbital and spin physics in LiNi02 and NaNi02 



3 



classical Goodenough-Kanamori-Anderson rules ^31 apparently predict to be FM. In 
a strong crystal field the ground state of Ni^+ is the low-spin [tlgc]^) configuration, so 
that direct exchange between t2g electrons does not occur. The exchange interactions 
between different nickel layers should also be FM ^1], so that one naively expects a long- 
range ordered FM ground state. However, for LiNi02 no long-range magnetic order was 
found down to temperatures very close to K and the magnetic susceptibility gradually 
diverges, giving the impression that the FM correlations mysteriously disappear. A 
suggestion was made by Feiner, Oles and Zaanen ^] that the eg orbital degeneracy of 
the Ni3+ ion and partly antiferromagnetic (AF) interactions might be responsible for 
this peculiar behaviour. 

The issue was then addressed by Mostovoy and Khomskii (MK) in an important 
paper j2] in which they proposed a realistic spin-and-orbital model for the Ni plane, 
which includes the Coulomb repulsion and the Hund's rule exchange splitting on oxygen. 
They arrived at the conclusion that there is a huge degeneracy in the orbital sector, 
which is not resolved at the mean-field (MF) level. Yet orbital order is favoured over 
an orbital liquid state by the order-out-of-disorder mechanism, while they claimed that 
anyway the magnetic interaction is always FM |2- From the absence of an orbital 
ordered state in LiNi02 they concluded that the difference between LiNi02 and NaNi02 
is probably extrinsic, due to disorder or electron-lattice interaction. Apparently this has 
now become the predominant view, and is as yet not inconsistent with experiments. 

However, a conclusion concerning the nature of the magnetic interactions and the 
origin of the peculiar properties of LiNi02 had better be drawn only after the theoretical 
prediction for the intrinsic in-plane behaviour is fully established. We believe that this 
is not the case and therefore reanalyze the situation in this paper. Our finding is that 
upon inclusion of the Hund's rule splitting also on the Ni ions and correction of what 
is apparently a mistake in the MK analysis, both FM and AF interactions can occur 
in the Ni plane, depending upon the orbital arrangement. Admittedly, this still leaves 
the difference between LiNi02 and NaNi02 to be explained, but it reopens the case for 
an intrinsic mechanism since different orbital phases in the Ni plane could give rise to 
different magnetic interactions. 

The paper is organized as follows. In section |2l we present the notation used in 
the following sections for the spin and orbital operators used to derive the microscopic 
model. In section 01 we review some issues concerning SE, in particular with regard to 
the Ni-O-Ni 90° bond, and compare this case with the standard situation encountered 
for 180° bonds in TM perovskites. The spin-orbital SE model for the triangular Ni 
planes of LiNi02 and NaNi02 is presented in section HJ Next we discuss the strong 
frustration of the orbital and spin interactions in this model, and we investigate its 
consequences and present possible ground states, obtained using MF theory for a pair 
of Ni ions (section Ej) and for the entire plane (section IH)). The implications of the model 
for the physical properties of LiNi02 and NaNi02 are discussed in section [7| Finally, 
in section |H1 the main conclusions and a summary are given. Technical details of the 



derivation of the model are presented in Appendix A 
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2. Pseudospin formalism for degenerate Cg orbitals 

For the twofold degenerate eg orbital state at each site local operators corresponding to 
pseudospin T = 1/2 are introduced, i.e. Tf, Tf and T/, acting as half the Pauli matrices 
(j^, cr^ and on the two-dimensional (2D) orbital Hilbert space at site i with basis 

(J) - V^) - 43.-.^|o), (°) - V-^) - <..-,.|o). (1) 

A general superposition is given by 

m = cos(0,/2) Q + sin(^,/2) (^^ , (2) 

which for example at 6*, = | corresponds to a d^^-xi orbital. The expectation values of 
the pseudospin operators in the general orbital state (j21) are 

(7;^) = lcos^„ (7;-) = isin^„ {Tr) = 0. (3) 

In order to make the formalism more flexible and to include explicitly the cubic 
symmetry of the Cg orbitals, it is convenient to deflne two more equivalent basis sets by 

lia) = dj,3a2-r2|0)' ({a,/5,7} acyclic 

\ia) = dl p2_y2\0) , permutation of {x,y,z}) (4) 

and corresponding rotated pseudospin operators and If behaving like Tf and Tf 
with respect to those basis sets, i.e. 

TX IT"'-^ \/3t^x tx _|__\/3t-iz 17^2; 

i 2 i 2 « ' i ' 2 i 2 i ' 

tV _ It^z _|_ VSrpx jy _ y/3^rpz iT^a; 

i 2 i ' 2 i ^ * 2 i 2 i ' 

I- = = Tr^, (5) 

which satisfy the identities 

i^ + iy + ^ = o, i^ + iy + ^ = o. (6) 

We can now introduce on-site orbital projection operators by 

pr = (ij, + jr), vf = {\i.-it). (7) 

where /j is the unit operator in the 2D orbital Hilbert space at site i. 

We further deflne two sets of (mutually dependent) orbital-pair operators, 
{X°^ Jtf} and {Vtf, Wtf}, all of which refer to a pair of nearest neighbour TM ions 
with their bond {ij) lying in the a(3 plane (see flgure^b)), 

= (1^ + vf){i, + ) + (/. + vf) (/, + Vf) 

= (f J. - /r)(|j, - If) + (§/. - If )(|/, - If), (8) 

= (/. - )(/, - vf) + (/. - vhiij - vf) 

= (ij, + imiij + /f ) + (ii. + /f)(ii, + If), (9) 

= -/.(/; + /f)-(/r + /f)i, = i./7 + /7/„ (10) 
w;^.'^ = 2 (If If + jf jj^) = 2 (J7 // -If If- If If). (11) 
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Their expectation values in a pair state \i0i)\j6j) are given by 

(X^f^) = |[35 + 12cos(^+ + Xf) cos^_ + 4cos2(^+ + x^) - 2cos2^_], (12) 

{J^f) = |[3 - 4cos(0+ + Xy) cose. + Acos\e+ + x-y) - 2cos2^_], (13) 

(Vf) = cos(^+ + x-y) cos^^_, (14) 

(>V°^) = i[4cos^(^^+ + X7) -2cos2^_ - 1] = i[2cos(2e+ + 2x-,) -cos2^^_], (15) 
where 9± = {6^ ± d,)/2, and Xx = ^iXy = -f ' X. = 0. 



Finally we introduce bond projection operators^ needed for specifying the SE 
interactions between a pair of TM ions. For the spin part we will use the familiar 
projection operators for spin triplet and spin singlet, 

Qij ~ 4 I'J + ' Qij ~ \'^ij ~ ' ^ji (^6) 

where lij is the unit operator in the four-dimensional (4D) spin Hilbert space on the 
bond {ij). For the orbital part we will make use of 

QZ% =prpf +pfp;, (17) 

Q-^^^. = V^V^ + Vfv^ + V^vf + pf pj^ , (18) 

QnI- = + Vfvf, (19) 

where the labelling will be explained in section |3] They are conveniently expanded in 
terms of the operators defined above, according to (with indices omitted for clarity) 

So = J = \I - \V + \-W 

= M - \Z - \J = I - W (20) 

Qn = -21 + iJ + = i/ + |V + iW , 

where / = 1^ now denotes the identity in the 4D orbital-pair Hilbert space. 



3. Superexchange for degenerate Cg orbitals 

If degenerate Cg orbitals are partly filled, a spin-orbital Hamiltonian can be constructed 
from a degenerate-band Hubbard model in much the same way as one derives the AF 
Heisenberg model from the single-band Hubbard model at half-filling. The Hamiltonian 
of the Cg-band Hubbard model consists of two parts: there is a hopping term T modelling 
transfer of electrons between nearest-neighbour transition metal sites and a Hubbard 
(or Coulomb) term U describing the on-site interactions. In the situation where the 
number of sites equals the number of electrons, the ground state for T = has twofold 
spin and twofold orbital degeneracy on each site. When we allow for a nonzero T as a 
small perturbation, the 4^-fold degeneracy is lifted by virtual electron hopping involving 
excited states, and the effective Hamiltonian in second order perturbation theory is 

Kff = E \r\n) ^ (n|r] Pg, (21) 

where Vq, is a projection operator onto the ground state manifold of W. The spin-orbital 
Hamiltonian is obtained by writing out ()21|) in terms of spin and orbital projection 
operators at each site i. 
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For a cubic lattice, such as in the perovskites KCuFa and K2CUF4, one can describe 
by this approach also the 180° SE, thus treating it formally as direct exchange. The 
hopping term is taken as 

T = -t E E E dL (22) 

Q (ij)\\a <r 

expressing that hopping is only allowed between nearest-neighbour directional orbitals 
\a) = c?3^2_r2|0) oriented along the connecting a axis (a being x, y or z). In the following 
we will therefore call the a-orbitals 'hopping' orbitals - since electrons in these orbitals 
can hop and contribute to the kinetic energy, and similarly call 'non-hopping' orbitals the 
a-orbitals {\a) = (i|;2_^2 10)), orthogonal to the a-orbital and oriented perpendicular to 
the bond. In a one- dimensional chain this situation leads to a characteristic competition 
between itinerant and localized phases |16j . 

The on-site interactions on a TM ion can be represented by [T7] 

WtM = Uj^^iX'i^iH + {U - -Ju) E '^iA'^iM - 2Jh J2 " 
iX i,\<fi i,X<fj, 

+ "^H E {dlx^dlxid'iij.idii^'^ + dl^^dli^^d^xidixi;^ J (23) 

i,X<fi 

and are characterized by two parameters: the intraorbital Coulomb energy U and 
the exchange energy Jh- The interorbital terms in equation (PHjl describe electron 
interactions between pairs of orthogonal orbitals, i.e., in the present subspace of eg 
orbitals one has A,/i G {a, a}. The excited states, generated in the virtual d — d 
transitions and relevant for SE, have two electrons on the same ion, and U contributes a 
Coulomb repulsion energy U. In the t <^ U limit one thus derives an effective low-energy 
spin-orbital Hamiltonian with coupling constant Jse i'^/U, in which spin and orbital 
degrees of freedom are interrelated. By taking also the Hund's rule exchange Jh into 
account one removes the classical degeneracy of magnetically ordered phases [TH]. The 
stable phase at low temperatures has long-range orbital order of a particular type of 
mixed orbitals, leading to AF interactions along the c-axis and FM interactions in the 
(a, 6)-plane. This ordering was verified experimentally and has been shown to be stable 
with respect to quantum fluctuations for large Jh [19j. 

In a more realistic treatment of SE one takes into account explicitly that the hopping 
takes place via the ligand oxygen ion. The hopping term is then taken as 

^ = E E E (4a<x Pjaa + p]aa diaa) , (24) 
« (ij)\\a ^ 

and describes charge transfer with amplitude t = t^j between the TM-orbital d^a'2_r2 
and the oxygen cr-type p-orbital Pa, where again the orbitals are oriented along the 
connecting a-axis. The on-site interaction on oxygen is given by 

Wo = f^O E "'jAT'^jAi + {Uo - ^ Jo) ^j^^jf^ ~ 2 Jo ■ 
jX j,X<fi j,^<fJ- 

+ ^ E {pUPmPj^^lPjf^^ + PmlP]f^lPj^lP,xi)^ (25) 
i,A</t 
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with intraorbital Coulomb energy Uq and exchange energy Jq. Here the operators UjXrj, 
Ujx and Sj\ refer to an oxygen ion at site j. For two- hole excitations in the three 
2p orbitals, as arises when the SE is derived (see below), this local problem and the 
excitation spectrum are isomorphic to those for three orbitals filled by two electrons 
as in the vanadates [20] • The effective Hamiltonian is now obtained in fourth order 
perturbation theory, 

Hefr= E Vo[T\k)-r^—r{k\T\l)-,^{l\T\m) ^ {m\r]Vo. (26) 

k,l,m^G ^G-H i^G-i^l ^G-^m 

In the case of a 180° TM-O-TM bond {ij) this is unproblematic. For the so- 
called [/-term (Anderson or delocalization process) J^l, where an electron is effectively 
transferred from one TM ion to the other TM ion, schematically represented as 

one can simply replace i by t^/A, where A is the excitation energy for transferring 
an electron from O to TM, and so the coupling constant in the effective Hamiltonian 
becomes 

*oc|i. (27) 

Note that the orbital filled in the second step is necessarily the same, also as regards 
spin, as the one emptied in the first step. The so-called A-term (Goodenough process 
or correlation effect) JH] involves instead electron transfer of two electrons from the 
connecting oxygen ion, one to each of the TM neighbours. 

Here the oxygen 2p^ configuration involved has two holes with opposite spin on the 
same a-type p-orbital, giving always the same intermediate state at the oxygen ion 
with energy 2 A -|- Uq, and the contributions to the effective Hamiltonian have coupling 
constant 

-^""^^^ ^2A + Uo ~ 2a) " 2A(2A + [/o)' ^^^^ 
The reason for the subtraction of the term oc 1/2A will be discussed below. For a 180° 
TM-O-TM bond the two processes (Anderson and Goodenough) make qualitatively 
similar contributions to the effective Hamiltonian (at least for a single eg electron on 
each TM ion j2I]) because both involve cr-type hopping. As the terms contributed to 
Tieff have identical form, they can be formally generated by the second-order formalism 
of direct exchange above, even though this models only the process giving the [/-term 
in the SE. To obtain the coupling constants one can simply add equations and 
(I28|) . As these are generally of the same order of magnitude, inclusion of the A-term 
is important quantitatively, and is essential to describe the trend in the strength of SE 
within the 3d TM series 
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In the case of a 90° TM-O-TM bond, as on a triangular lattice, the situation is 
very different. In the [/-term (Anderson) process the oxygen orbital through which the 
electron is being transferred is now cr-type for one TM ion but 7r-type for the other one. 
Therefore this process can only occur if it involves a (deep-lying) t2g orbital, 

{4e')p%ty) - {4e')p'{4e') - {tle')p%ty) , 

and thus contributes to the effective Hamiltonian a term with coupling constant 
t'^t'^ 1 

since the energy of the middle intermediate state is increased by the crystal field splitting 
AcF, and the hopping parameter for the t2g orbital is instead of t^r- In the charge 
transfer terms (A process) the oxygen 2p'^ states now involve two holes on different p- 
orbitals, each of a-type but with respect to a different TM neighbour. This implies that 
the oxygen Coulomb interaction involved is now the interorbital interaction Uq = Up 
instead of the intraorbital interaction Uq = Uq + 2Jp. Moreover it leads to a singlet- 
triplet splitting, with energies Us = Up + Jp and Ut = Up ~ Jp, where Jp = Jq 
is the (Hund's rule) exchange on oxygen, and thus the contributions to the effective 
Hamiltonian have coupling constants 

j± ( ^ M - Up^Jp 

^ A2V2A + f/p±Jp 2A^ A2 2A(2A + ?7p± Jp)' ^ ' 

From equations (j29|) and (j3(J|) one observes that the dominant contribution to 
the SE comes from the A-term (Goodenough process), since t^/t^ ~ 4 ^24], while 
U + AcF ^ 2 A and Up ~ 2 A [23 1211 l2Zj • This was already pointed out in reference 
[13] , as well as the fact that the magnetic interaction should then be FM, as equation (jHn|) 
shows, which is one of the famous Goodenough-Kanamori- Anderson rules. Nevertheless, 
one should be careful not to jump to conclusions here, since the implications of the 
orbital SE interaction were not fully considered by Goodenough (the case explicitly 
considered was the SE between Ni2+ (t^^e^ : M2) ions [2H1, where the two Cg orbitals 
are both occupied by one electron, and no orbital effects can arise). It is further clear 
that the dominant A-term (Goodenough contribution) cannot be represented well as an 
effective second order direct exchange. Yet this was attempted in a recent paper P^j : 
this has the merit of deriving the most general form of the effective Hamiltonian purely 
based on symmetry considerations, but it cannot capture the dependence on the most 
relevant parameters Up and Jp (actually only the weaker U-term SE, dependent upon 
the splitting of the intermediate Ni^"*" configurations, is being described). 

A notable peculiarity of the A-term, already included in equations ([2H1) and ([301), 
is that subtraction is needed of the contribution that would have been obtained if the 
electrons transferred to the TM ions would have come from two different oxygen ions 
and not from the connecting oxygen ligand, as pointed out by Mostovoy and Khomskii 
[21]. The necessity for this subtraction can be understood as follows. The reference 
state ('vacuum') should be considered to be renormahzed by all possible fourth order 
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Figure 2. Double charge transfer excitations in LiNi02: (a) Isolated excitations which 
renormalize the electronic structure, and (b) Joined excitations involving a common 
oxygen site, which contribute to the superexchange. 

uncorrelated hopping sequences [see figure Efa)], each contributing a term t^/(2A^). In 
most cases they cancel out because the three contributions corresponding to hops 2 and 
3 in figure I2^a) made along the three cubic axes with the axis of hops 1 and 4 kept 
fixed, add up to a constant, since the three projection operators do so because of 
equation ©, and thus only add to the vacuum energy. However, the cancellation fails if 
the contribution from one axis is missing because the oxygen ion there is joined with the 
other TM ion, as in figure |2l^b). The cancellation is restored by adding this term to the 
other two and subtracting it from the SE for the pair of TM ions under consideration. 

Note that the above correction implies a sign change of the coupling constant and 
so the opposite situation is favoured than one might naively expect. In particular, 
the largest diagonal SE is generally obtained for the configuration that permits the 
largest number of hopping sequences returning to itself, and so this configuration is 
now energetically disfavoured instead of favoured. This applies specifically to purely 
interorbital SE terms, whereas in spin-spin SE interactions, which generally originate 
from the difference of SE for two spin multiplets, the corrections cancel out and one 
gets the expected result. 

4. Spin-orbital model for LiNi02 and NaNi02 

Let us now reconsider the derivation of the spin-orbital orbital model for the triangular 
lattice structure of LiNi02 ISO], taking only the charge transfer process (A-term) into 
account, as done also by MK So we consider two nearest neighbour Ni^"*" ions 
connected by two Ni-O-Ni 90° bridges, as in figure [H or |21^b), and analyze the fourth 
order hopping sequences (which in this order of perturbation theory can be done for 
each bridge separately). In order that these virtual excitations lift the ground state 
degeneracy it is essential that the intermediate states are affected by the on-site Coulomb 
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Figure 3. Possible charge excitations from the e^p^ to the e^p^ configuration. The 
oval represents the TM ion with one Cg electron, either in the non-hopping orbital 
d^2_^2 on the left or the hopping orbital d3a2_j.2 on the right. The circle represents 
the oxygen ion with the electrons in the orbital pa in the segment on the left. The 
excited states '142, and ^Ai at the TM ion reached in each case by the u-type 
hopping are labelled by T, D and S, respectively. 

interactions on oxygen and/or nickel, described by U, see equations and (j2Sl)- 

As indicated above, the relevant states of the oxygen p^ configuration, i.e. those 
occurring upon hopping, are a triplet ^Ti and a singlet ^T2, denoted for brevity by t 
and s. They are split by 2Jp, while moreover the interorbital Coulomb repulsion Up for 
p^, absent for two p^ configurations, must be taken into account. The relevant Ni 
states are Mi, and for which we will use the abbreviations S (singlet), D (orbital 
doublet) and T (triplet), respectively. These terms have equidistant energy levels given 
hj U + Ju, U — Jii and U — 3 Jh, with the triplet being lowest by Hund's rule, where the 
Ni^"*" Coulomb repulsion and Hund's rule couphng can be expressed in terms of Racah 
parameters asU = A + AB + 3C and Jt, = AB + C ^Hj. 

The first transition in the hopping sequence is e^p^e^ e^p^e^, with excitation 
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Figure 4. Example of a double charge transfer excitation e^p^e^ e^p'^e^ involved 
in the superexchange. The triplet and singlet states at oxygen are labelled by t and s; 
the labels for the states at TM ions are the same as in figure |3| 

energy 

Ax = Ef{e^) - E^Xe') + - ^°(/), (31) 

i.e. the charge transfer energy depends upon the Ni^+ state accessed, 

As = A + 2Jh, Ad = A, At = A-2Jh, (32) 

where it is understood that U — Ju has been absorbed into A. Which states are reached 
depends on the Cg electron aheady present on the Ni^+ ion to which the electron hops, 
as illustrated in figure which shows the six possible hopping channels. If the d^a'^^r'^ 
orbital is occupied, then a hop from the Pa oxygen orbital is only possible for an electron 
with opposite spin, and the excited states involved are the spin singlets Mi and ^E. On 
the other hand, if the df^2_^2 orbital is occupied (and therefore the d3a^_j.2 orbital is 
empty), then the spin of the hopping electron can have either sign, and the and 
states are reached. In the second step e'^p^e^ e^p'^e^ , the excitation energy is raised 
further by Ay + f/p ± Jp, depending upon the state Y accessed at the other Ni ion 
and the p^ state (s or t) at the oxygen. Figure El shows an example of such a double 
charge transfer excitation. In the third and fourth step the excitation is undone, either 
in the same or in reverse order. 

To derive the complete spin-orbital Hamiltonian one must list all possible initial 
configurations, and for each of them list all possible hopping sequences that return to the 
ground state manifold. The initial and final states are described by means of projection 
operators, both for the orbital occupation and for the spin state, defined in section |21 
The orbital bond projection operators specify whether the 'hopping' orbitals on the two 
Ni^"*" ions are both occupied by an electron (a situation denoted by 'O'), whether one 
hopping orbital is occupied while the electron on the other ion is in the non-hopping 
orbital ('M', for mixed), or whether both non-hopping orbitals are occupied by the two 
electrons ('N'). Since the contributions from the two Ni-O-Ni bridges are independent 
and may be simply added, the operators (fTTHT^ are defined to do so for each a/3 bond 
direction. Obviously these operators depend on the bond direction, and so this gives 
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explicit orbital anisotropy. For specifying the overall spin state (i.e. formed by the two 
spins 1/2 of the two Ni'^+ ions), we need the familiar projection operators ()16p for spin 
triplet and spin singlet. 

Collecting the contributions from all configurations and all hopping sequences with 
the help of diagrams like in figures El and |3 one then obtains the Hamiltonian in first 
instance in the form 

KfT = E {Qt,lKoQl+Kl>Ql]+QtMQl+K^^^^ (33) 

where it is understood that the form of the projection operators depends on the bond 
direction. In order to separate the spin dependent part from the purely orbital part this 
may be rewritten as 

'^eff = E ( ["^O + "^M QM,ij + "^N Qn%] '^ij + ["^O + QM,ij + Qn%] '^j)^ (34) 

with orbital and spin-orbital interactions 

Jl = ^,Kl + lKl Jl = Kl-Kl (L = 0,M,N). (35) 

To see how this works in detail and to compare with the analysis of MK, let us 
initially ignore the Hund's exchange splitting on Ni (the full procedure is explained in 
more detail in Appendix A). If the initial configuration is N-type, the electrons hopping 



from the different 2p orbitals at the common oxygen ion into the empty hopping orbitals 
on the Ni neighbouring ions, can do so with their spins oriented in four ways: both up or 
both down, thus leaving the oxygen in the triplet p^ state, or with their spins up-down 
or down-up, both with equal probability 1/2 for leaving the oxygen in the triplet or 
in the singlet state, all in all making three triplet and one singlet contribution. As 
the terms on the Ni ions are all equivalent when Jh = 0, this is independent of the 
initial orientation of the spins in the Ni non-hopping orbitals. Upon including an overall 
factor of 4, because both the two excitation transfers and the two deexcitation transfers 
can also be made in reverse order, one obtains 

Kl = Kl^ = l2[XtX] + A[XsX]. (36) 

Here we denote the fourth-order perturbation expressions by giving a shorthand notation 
for the middle intermediate state (using here X instead of D or T, because we do not 
distinguish between those states yet). Their values are [compare equation (jHn|) above] 

[^^^] = ( " 2A + f/p - Jp ^ 2a) " 2A(2A + f/J- J,V ^^'^^ 



p) 



\X^X] = — ( ^ I M = Up + Jp /ooN 

^ ^ J A2 V 2A + t/p + Jp 2A^ A2 2A(2A + t/p + Jp) ' ^ ' 

If the initial configuration is M-type, the result is equally independent of the initial 
Ni-Ni spin state. Although the spin of one transferred electron must be opposite to the 
spin of the electron occupying the hopping orbital on Ni, the other transferred electron 
can have its spin still either parallel to that of the first one, leaving a triplet state 
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on oxygen, or antiparallel, yielding a triplet or a singlet with probability 1/2 each. 
Including again the factor 4, one obtains 

K^ = Ki, = 6 [XtX] + 2 [XsX] . (39) 

Only if the initial configuration is 0-type, the result is different, because each of 
the transferred electrons must have its spin opposite to that of the electron already 
occupying the hopping orbital on Ni. So, if the spins of the electrons on Ni are parallel, 
the initial state thus being a Ni-Ni spin triplet, then the transferred electrons necessarily 
also have parallel spins, leaving oxygen in the triplet p^ state. If the electrons on the 
Ni ions have antiparallel spins, i.e. being either in a spin triplet or in a spin singlet 
state depending upon the phasing between up-down and down-up, then the spins of the 
electrons being transferred are also antiparallel, and have the same phasing, so again 
the Ni-Ni spin triplet yields an oxygen triplet, while the Ni-Ni spin singlet yields an 
oxygen p^ singlet. Therefore, upon inclusion of the factor 4, 

= 4[XtX], = A[XsX]. (40) 

It follows that 

J°=4Jr, A = 2Jt, J'o = Jt, (41) 

4 = 0, 4 = 0, JE = -Jts, (42) 

J - 3\XtX]-,\XsX] ^(277, - J,) + 7/,^ - J,^ ^ 2t^ U, 

Jt - 3[XtX] + [XsX] -— ^2A + f/,)2-j2 -A^^ATf^' 

J.S - -mm - [XsX]) = % ^ % j^^2, (44) 

where the final expressions on the right in equations (pSj) and (jl^ are the results in the 
limit Jp <^ A, Up. So the effective Hamiltonian is 

{ij) 



{ij) 



"^lij — Jts Jif Si ■ Sj], (45) 



with 2°'f and J'^^ given by (jHl) and Q, and understood to depend upon the direction 
of the bond {ij). 

Comparison with the results of MK |2j shows that our analysis confirms their finding 
that the purely orbital interaction is stronger by one order of magnitude than any spin 
dependent interaction, ^ Jts, in agreement with the conjecture made by Reynaud 
et al and based on their experimental data JU]- One recognizes that this comes about 
because, at Jh = 0, all spin dependence originates from the singlet-triplet splitting of the 
oxygen 2p^ configuration, and Jts is therefore smaller by a factor ~ Jp/Up ^0.1 than the 
orbital interaction Jt. Our analysis also confirms the form of the orbital interaction as 
being oc l^f ■ However, we find a different form for the orbital dependence of the mixed 
spin-orbital interaction, since MK give X°-^ instead of Jtf ■ Apparently their result 
is incorrect, because the reasoning above clearly demonstrates that only the 0-type 
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configuration can give spin dependence, and this is what is expressed by the operator 
j-ai3 _ -p^^ The difference is important, because the bounds on the expectation values 
are very different, y > (T"'^) > y whereas | > {J"^) > [see equations (fT^ and 
(fT^]. and so MK concluded that the spin exchange is effectively always FM for any 
orbital state, while we conclude that this is at best marginally so, i.e. the spin exchange 
interaction on a bond {ij) given by ()45|) can vanish for some specific combination of 
orbital states at sites i and j. 

Moreover, the latter conclusion is only valid up to the present level of 
approximation, where only the Coulomb interactions on oxygen have been taken into 
account. If we include also the Coulomb interaction on nickel, i.e. allow Jh to be finite, 
and consider the leading correction to (03), we find that the spin exchange may even 
become weakly AF. To see this we first observe, as demonstrated in [Appendix A that 
quite generally the coupling constants satisfy 

I J^l > Ji^ > I Jj^l withJi^>0, (46) 

so that the next spin-dependent term to be considered is 

{ij) {ij) 

and next that it follows from equation p5|l that | > (/ — W"'^) > |. We will investigate 
below under which conditions an AF net spin exchange on a bond (ij) could actually 
occur. 

For that purpose it is useful, when allowing both Jy^ and J§ to be finite, to rewrite 
the full spin-orbital SE Hamiltonian (j3l|l explicitly as the sum of a purely orbital part 
and a spin dependent part, 

HcS = 'HcS,o + '^^eff.s, (48) 



Hefr,o = E [J'tW + JtJ^' ) 1^.- = JrY: y^tf 1.., (49) 

{i,j) {iJ) 

^cfr,s = E [j'TsIij - J'TS^if - JTsJtf ) Si ■ Sj, 



{id) 



E ( - J<rh, + JuVtf - J.Wtf ) 5, ■ 5„ (50) 

{i,j) 



with the exchange constants all positive and satisfying (see [Appendix A ) 



< < < < Jr. (51) 

Note that is the pure spin SE, while Jy is the strength of the SE between a spin-and- 
orbital operator I'^Si on one site and a pure spin operator IjSj on the other site and 
is the SE between spin-and-orbital operators on both sites. So equation shows that 
the pure spin SE is actually always the weakest interaction (for the above case of Jh = 0, 

one has = Jt, Jt = 0, Jr = \Jt\ J'ts ~ J'ts ~ 0' ^ts = Jtsi Ja = J^i = Ju = \Jts)- 
In the first expression in equation ()49|) we have dropped a constant, while in obtaining 

^(«.i> ^ij 



the second expression we have used that Y^u j) '^tf = 0; follows from equations (fTUI) 
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Figure 5. Interaction constants (in units of J = t^/A^) as a function of Up/ A for 
Jp/Up — 0.1 and Jh/A — 0.05: (a) Orbital interaction constant Jr] (b) Spin-orbital 
interaction constants J^, Ji, and J^- (from top to bottom). 

and (jni). Note that the latter simphfication is not possible in equation dSDI), where the 
spin inner product occurs instead of the unit operator in spin space. 

To get an idea of the order of magnitude of the SE coupling constants one would 
need to estimate the parameters involved. Unfortunately, reliable estimates are available 
almost exclusively for the divalent TM ions, while only little work has been done on the 
trivalent ions fIF)\ EEl l2Zj- Based upon what is available we estimate (all values in 
eV): t ~ 1.5, Jh ^ 1.1, Up ~ 5, Jp ^ 0.8 [S!] and A ~ 2, where the last value is 
the most uncertain. The small value of A suggests that the covalency effects are quite 
important in LiNi02, and thus the present model has to be considered only as describing 
the generic structure of the interactions between states with spin S = 1/2 and e^-type 
orbital degeneracy, which originate from Ni^"*" ions surrounded by a hole shared between 
the Ni 3d orbitals and the oxygen 2p orbitals. It is anyway obvious that A is too small 
to consider the perturbation theory a controlled expansion. We shall therefore regard 
the overall energy scale J = t^/A^ as an unknown parameter and take it as our energy 
unit. We can then study the relative size of the coupling constants, which are controlled 
by the ratios Up/ A, Jh/A and Jp/Up, which we treat as parameters. Where any of 
them needed to be fixed for varying other variables, we took Up/ A = 2, J^/A = 0.05 
and Jp/Up = 0.1 (which is maybe a bit small for the first and especially the second and 
third one). 

The variation of the orbital SE coupling constant Jr and the three spin and spin- 
orbital SE constants J^^, J^ and J^ with the strength of the oxygen Coulomb repulsion Up 
is shown in figured It illustrates once again that oxygen Coulomb interaction is crucial: 
at Up = all coupling constants vanish. One further notes that the inequality (jHT|l is 
well satisfied: indeed Jr is by far the largest, and the three spin coupling constants are 
in the expected order. Note also that there is more or less saturation for Up/ A larger 
than ~ 2 so that the inaccuracy in this parameter is not so important. 

Figure El shows the variation with Jh, the strength of the Hund's rule exchange on 
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Figure 6. Interaction constants (in units of J = i^/A^^) as a function of 2Jh/A for 
Up/A = 2.0 and Jp/Up = 0.1: (a) Orbital interaction constant Jt-; (b) Spin-orbital 
interaction constants J^^, Ji, and Ja (from top to bottom for positive Jh)- 

nickel. One observes that a finite Jh enhances the orbital SE considerably [figure IHl^a)], 
basically because this lowers the energy of the triplet and so further disfavours (see 
the remark at the end of section ^ the N configuration, already disfavoured because it 
allows the largest number of hopping sequences. Moreover, figure Efb) shows explicitly 
that nonzero Jh is essential for having unequal spin interaction constants. Since this 
is equivalent to Jj^ being nonzero (and positive) and thus for having potentially AF 
interaction, as argued above when the Hamiltonian (jTfjl was discussed, inclusion of the 
Ni^"*" term splitting makes a qualitative difference for the obtained spin-orbital model. 
It is amusing that in the unphysical regime of negative Jh the order of the spin SE 
constants is inverted, so that J^ would be negative and AF interactions would not be 
possible. Note further that magnetic frustration is apparently maximal at Jh = 0, where 
the three types of spin dependent interaction are equally strong and so experience the 
strongest mutual competition. 

5. Interaction between a pair of Ni^+ ions 

In view of the dominating size of the orbital interaction J^- it is natural to consider first 
the orbital interactions alone. Before addressing the Hamiltonian ()49|1 on the entire 
triangular plane, it is useful to solve first the problem of a single Ni^+-Ni'^"'" pair, which 
for definiteness we assume to have its bond in the XY plane. For an arbitrary orbital 
pair state \i6i)\j6j) the orbital energy of the pair is then given by [compare equation 

Eri0i,9,) = ij,[4cos2(^) - 2cos2(^) - 1], (52) 

and the corresponding energy surface is shown in figure [7|^a). One notes that similar 
orbit als {6i ~ 6j) are generally favoured over dissimilar orbit als {6i ~ 6j ± tt). More in 
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Figure 7. Orbital superexchange Jt for a single pair along XY (see figure^: (a) 
bond energy as a function of the orbital angles 6i and 62 (energy in units of jJt, angles 
in units of tt); (b) the two favoured pairs of occupied Cg orbitals. 



particular, for a single bond 7^efr,o is seen to favour specific identical orbitals, viz. those 
with 9i = 9j = i.e. a pair configuration with orbitals {\z)±\z)) / \^ oc d3z2-r2±dx2-y2 
at both sites, as illustrated in figure Ulh), with energy Ej.^" = — |Jr- Frustration 
prevents this optimum orbital arrangement to be realized on all bonds simultaneously, 
as already pointed out by MK. As will be seen in the next section, on the full triangular 
lattice the orbital interaction favours phases in which a pair has either a ferro orbital 
(FO) arrangement at a general angle 6 = 6i = 6j or similarly a canted orbital (CO) 
arrangement with 6 = 6i = —6j. It is therefore of interest to consider the effective 
spin-spin interactions realized under those conditions. 

Quite generally, once the orbitals of the pair are fixed, supposedly determined by 
the pure orbital interaction, the resulting effective spin SE interaction can be written 

Heisenberg Hamiltonian, 



cff/ 
ij \ 



i,Oj 



with the effective spin SE coupling J^f given by 



Si ■ Sj, 



(53) 



j: 



eff/ 



^J y^^,^JJ ^ ^4 - (l + 2 COS^ 9^-4: COs'^{9+ + Xy)) - "^Jv COs{9+ + Xy) COS 9_] , (54) 

where 9± = {9i ± 6'j)/2, as obtained by inserting equations (fT^ and (fT3j) into (jHOJ. For 
the two relevant cases of orbital order one thus finds, by adopting the appropriate values 
for the angles 9i and 9j, 

(55) 
(56) 



JcoiO) = J\ 



eff/ 
ij \ 

cff/ 



4J^ - 3 + 4J^ cos^(^ + Xi) - "^Jv cos{9 + ^7)] , 
I 4Jo- — (1 — 4 cos^ X'y)JtJ. ~ 2 cos^ 9 — 4: cos x-y^i^ COS 9 

The dependence on 9 is shown in figure |H1 for the values of Jo-, and obtained 
with our standard parameter set. Figure Ufa) shows that in the FO case the spin 
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Figure 8. Effective spin-spin superexchange constant (in units of J) versus angle 
6 (in units of tt; see equation Q): (a) For a single pair along XY for ferro orbitals 
(0 — 9i — 9j); (b) For a single pair along YZ or ZX for canted orbitals {9 — 9j — —Oj). 

SE constant is mostly FM (negative), as expected from the Goodenough-Kanamori- 
Anderson rules, but that AF values, though smaller, are indeed possible, notably for 
orbitals with 6 ^ ±7r/3. Note that the coupling constant is that for a bond along XY; 
equation shows that the curve should be shifted by ±27r/3 for bonds along YZ or 
ZX. The CO case is shown in figure IHl^b), for a bond in the YZ or ZX direction while 
the angle is still given with respect to the basis {\iz), \iz)} as in equations (P) and 
Note that the coupling can again be AF as well, precisely for those angles {9 ^ ±7r) 
where the FO arrangement gives the largest FM coupling in the XY direction. 

6. Possible ground states 

We have already noticed that the orbital SE is strongly frustrated — since the 
specific orbitals minimizing the pair energy are different for each bond direction, the 
interactions between neighboring Ni ions cannot be satisfied in all three bond directions 
simultaneously. The fact that frustration occurs although the interaction is of ferro 
type is somewhat unusual and is generic for orbital physics as compared to spin physics. 
The underlying reason is that the orbital interaction Hamiltonian pUj) does not have 
global continuous rotational U{1) symmetry with regard to the orbital angles 6i, but is 
only invariant under rotation by an angle ±27r/3, corresponding to a cyclic permutation 
of the bond directions XY, YZ and ZX. This characteristic feature and the ensuing 
frustration are not captured by S'f/(4) symmetric models |S21E3- 

Next we consider what states with long-range order (LRO) are possible. As pointed 
out by MK, the simplest such state, viz. with uniform ferro orbital order, is actually a 
minimum energy state in MF theory. Its energy is independent of 6 and equal to — | Jr 
per site, while moreover this FO phase may be interrupted without energy cost by lines 
in which all orbitals have the opposite angle —6, i.e. by canting of all orbitals along a 
line in the triangular lattice. In order to see how this comes about, one may write down 
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the energy for a set of FO lines (for definiteness along the XY direction) with orbitals 
in line n given by On- From and (fT3j) one obtains 

-^FO— 1 -^intraline ~l~ -^interline 

= iJ,5][2cos20„-l] 

n 

+ lJr E [cos {On + On+1 " f ) + COS {On + ^„+l + f ) " COS {6^ ~ 

n 

= l-JrYl[4: COS^ 9n - 3] - JrY^ COS 6'„ COS 6'„+i (57a) 

n n 

= iJrE [2(cosft„ -cos0„+i)' -3]. (576) 

n 

Equation (j57aj) shows that upon variation of the angles On the gain in intraline line 
energy is compensated by a loss in interline energy, leading on balance to Epo-i being 
not dependent on the value of 9 in the ordered phase, while equation ()576|) shows that for 
any two neighbouring lines the same energy is obtained for 9n = +dn+i and 6'„ = —9n+i- 
So in addition to the degeneracy with respect to 9 generated by the frustration, the 
model also gives rise to randomness due to uncorrelated switches from +9 to —9 between 
successive lines. Such a random CO phase is probably best interpreted as the occurrence 
of twin boundaries or antiphase boundaries between equivalent FO domains. Because 
of the zero energy cost for formation of a boundary there is actually no preference for 
the formation of large domains. 

It was argued by MK that the degeneracy in 9 is an artefact of MF theory, and 
that this U{1) symmetry of the FO state, not present in the original Hamiltonian ()49|). 
would be lifted in the FO domains by the 'order-out-of-disorder' mechanism [H^, and 
that in this way the threefold symmetry associated with the triangular lattice would be 
restored. However, as the energies involved in the restoring quantum fluctuations are 
generally small, only a fraction of the orbital coupling constant J^, one should at this 
stage reconsider the orbital-spin coupling. Even though the relevant coupling constants 
Jcr, Jfj. and Jj, are much smaller than Jr, the spin-and-orbital dependent SE energy might 
still be more significant than the fluctuation energies in lifting the 9 degeneracy. Figure 
IE] suggests that it could do so already at the MF level and so determine the ground state 
before fluctuations need to be considered. 

We therefore consider the following LRO spin patterns in conjunction with orbital 
order: (i) FM order; (ii) the 120° three-sublattice order which is the ground state of the 
classical TALAF; (iii) AF line order similar to that of the orbitals, i.e. lines of parallel 
spins along the same direction as the orbital lines, and the spins alternating between up 
and down on successive lines. The result for the MF energy is shown in figure Efa). It 
turns out that indeed the 9 degeneracy is lifted, and the lowest MF energy is seen to be 
obtained for FO order with orbital angle 9 = +7r (or 9 = —tt), corresponding to a d^i-y^ 
orbital at every site, in combination with spin AF line order, illustrated in figure IHl^b). 
This shows that even though the spin-orbital interaction is considerably weaker than 
the pure orbital interaction, it can still have a qualitative effect on the nature of the 
ground state by breaking the U{1) symmetry generated by the orbital interactions alone 
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Figure 9. Spin-and-orbital order: (a) MF energy per site (in units of J = t'^/A^) 
for various ordered phases as function of the orbital angle 6: (upper horizontal line:) 
ferro-orbital / 120° spin order ; (lower horizontal line:) ferro-orbital / ferromagnetic 
order, degenerate with orbital-line / ferromagnetic order; (curve:) ferro-orbital/ spin- 
line order, degenerate with orbital-line / spin-line order; (b) Ground state in mean field: 
ferro-orbital order of (1^2 _y2 orbitals with spin- line order (along the XY direction). 

and restoring threefold symmetry. Therefore the large difference in coupling strength 
does not necessarily lead to a decoupling of the spin and orbital degrees of freedom. 

Remarkably, the +9 / —9 twin-boundary degeneracy persists, not only at the 
optimum angle ±7r, where it is immaterial since it corresponds only to a sign change 
of the orbital wavefunctions, but also at arbitrary 9. This could in fact have been 
recognized from the 9 dependence of the effective spin exchange constants: the spin 
exchange in the transverse direction (on the bonds along Y Z and ZX) as given for the 
FO case by shifts over ib^ in figure IHl^a) is the same as that given for the CO case in 
figure IHfb). 

Most significantly, we note that both for the ordered FO phase and for the 
disordered CO phase the effective spin-spin exchange is antiferromagnetic in the 
transverse direction, also for orbital angles deviating somewhat from the optimum 
value dzTT. However, for a disordered phase with truly random orbital angles, either 
a paraorbital state at high temperature or an orbital glass, the situation would be very 
different. Figure |H1 indicates that orbital randomness would produce a wide distribution 
of mostly FM exchange constants, although small AF values would still occur but not 
in a regular pattern as in the FO and CO phases. Obviously a phase transition to an 
orbital ordered phase would manifest itself also in the magnetic properties, e.g. in the 
dispersion of spin waves. 

Whereas the above analysis shows that the spin-orbital interaction could support 
orbital LRO, the degeneracy in the orbital sector and the small MF energy of only — | Jt 
per bond, both for the FO and for the CO phase, raise the question whether orbital LRO 
by itself is stable against formation of a disorded state of valence bond (VB) type. We 



Orbital and spin physics in LiNi02 and NaNi02 



21 



3. 




Figure 10. Geometry relevant for the orbital valence bond (OVB) state: (a) 
Triangular three-site plaquette in a Ni plane for determination of the exact eigenstates 
of the orbital interaction; (b) Lattice covering with thrcc-sitc plaquettes in the OVB 
state. 

therefore studied orbital valence bond (OVB) states [1^1 E2] in the following way. First 
we solved the eigenvalue problem of the orbital SE Hamiltonian 7ieff,o exactly on the 
three-site triangular plaquette shown in figure ITUT a). The eigenstates are four doublets, 
with ground state energy i?Q^^^ = — and excited states at — Jr, +^Jt and +2Jr, the 
ground state doublet being 

= ^{\x)i\z)2\y)3+ \z)i\y)2\x)3), 

\e2) = ^{\x)^\z)2\yh + \z)i\y)2\xh). (58) 

Note that particular combinations of orbitals are favoured, similar to the orbital order 
for the V triangle in LiV02 [33]. It is remarkable that on a triangular lattice such 
triangular VB correlations in orbital space are favoured, in contrast to dimer singlet 
correlations in spin systems jSHj- 

We then constructed an OVB solid by covering the triangular lattice with three-site 
plaquettes, as in figure ITUTb). and assigning to each triangular plaquette / an orbital 
wavefunction built from the groundstate components (j58|l . 

|A,) = cos6', |ei)i + sin6'i |e2)i, (59) 

with 6i still to be determined. This amounts to assigning a fictitious spin to each 
plaquette and solving a frustrated spin problem on the triangular lattice. In MF theory 
the total energy is now given by addition of the groundstate energies from all plaquettes 
and the interplaquette energies coming from the two bonds between nearest neighbour 
plaquettes, 

EoYB = E ^A, + E(^^ A^|Hefr,o| Ai, A,„) = E ^o'"" + E Ey'B'idl^ ^rr.). (60) 
/ l,m I l,m 

For a uniform solution the energy per bond is given by 

EovB = 1(2 X I X ^o'^' + 4 X i X E^g') , (61) 
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by counting the number of intraplaquette and interplaquette bonds. Although the 
plaquette eigenenergy amounts to ^Eq^^'^ = —\Jt per bond of the plaquette, which 
for this exact solution is of course lower than the energy per bond for the FO and 
CO phases, the plaquettes cover only one third of the bonds of the lattice, and so the 
contribution from the interplaquette energy is decisive. 

The optimum values of the plaquette angles are found to be 6^/ = 6'^ = — ^, for 
which E^~^ = —^Jt, and if this value could be attained between any pair of plaquettes 
one would have -Eqvb = — (| + -^)Jt — —0.35037,-, lower than Epo = Eqq = — 0.25J^. 
However, this would requiring alternating angles which is impossible because of 
frustration on the triangular lattice. The optimum value that can be obtained is 
E'ovB — — 0.228Jr, realized for instance by line order of ±|. This is remarkably close to 
the energy per bond of the site-ordered states, and just confirms that the orbital sector 
is strongly frustrated. It is therefore likely that some quenched disorder would suffice 
to turn the system into an orbital glass below a freezing temperature T^i determined by 
J^, as proposed to explain the behaviour of the orbitals observed in LiNi02 [TUj . 

7. Origin of the difference between LiNi02 and NaNi02 

In view of the above results the properties of LiNi02 are puzzling CHI, as neither 
orbital nor spin order sets in down to the lowest temperatures. Thus the high degeneracy 
of the ground state manifold apparently persists, with Ni^+ ions being in the low-spin 
(^2g^g) state, with twofold orbital and twofold spin degeneracy. The standard scenario 
for such ions would be a cooperative Jahn- Teller effect lifting the orbital degeneracy 
below a structural transition at T^, followed by a magnetic transition at T/v, lifting also 
the degeneracy in spin space. Both transitions could be induced already by the electronic 
mechanism involved in the SE, but the interactions with the lattice are expected to play 
a significant role in the structural transition, enhancing the value of T^, as found in 
the manganites ||22|- In spite of the strong orbital interactions reported above, and 
although the S = 1/2 spins interact by the predominantly FM SE described in the 
previous sections, and the magnetic susceptibility follows the Curie- Weiss behaviour, 
X oc (T — 6'cw)~^ with 6'cw ~ 35 K down to about 80 K IH El E], long-range order 
in either orbital or spin space is absent [311]. This would suggest that the Cg orbital 
degree of freedom plays an important role, possibly stabilizing strong singlet orbital 
correlations on individual bonds, and supporting a spin liquid state [HI . 

Before addressing the central question of the origin of the spin liquid state in LiNi02, 
let us summarize shortly the structural and magnetic properties of the similar NaNi02 
compound, which crystallizes in the same structure as LiNi02 (see figure The 
standard scenario with two subsequent transitions described above is indeed realized in 
NaNi02. The latter compound undergoes a first-order cooperative Jahn- Teller transition 
lowering the local symmetry from trigonal to monoclinic at Tg = 480 K jHE], and the 
Ni-0 distances change from d = 1.98 A to two long bonds with d = 2.14 A and four 
short bonds with d = 1.91 A at T < Tg. In the low-temperature phase a magnetic 
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transition at T/v = 20 K to an A-tjpe AF insulator follows — it was first derived from 
magnetization measurements on a single crystal by Bongers and Enz already in 1966 
|39j ■ This magnetic phase has also been confirmed by a complete static and dynamic 
magnetic study at low temperatures, where strong anisotropy between the effective 
magnetic interactions was established jHE]- An FM intralayer coupling of Jfm = —26 
K and a weak AF interlayer coupling of Jaf = 2 K were found, values consistent with 
the observed Curie- Weiss temperature 9 = 35 K and Neel temperature = 20 K. 
Similar values of the exchange constants (Jfm — —29 K and Jaf — 1-9 K) were also 
deduced recently by analyzing the dispersion of spin waves found in neutron scattering 
experiments jlOj. We remark that the parameters which fit the Curie- Weiss law change 
around the structural transition at Tg fSH'], which indicates that the magnetic couplings 
depend on the actual orbital state, in agreement with the model presented here. 

In contrast, in LiNi02 the structural transition is absent, but EXAFS experiments 
indicate the presence of local Jahn- Teller distortions below an orbital freezing 
temperature Tof ~ 400 K, with two different Ni-0 distances again favouring occupied 
directional d^z^_j.2-tYpe orbitals [H]: two long bonds with d = 2.09 A and four short 
bonds with d = 1.91 A. While these local distortions are similar to those observed in 
NaNi02, the absence of a macroscopic distortion in LiNi02 has been a mystery since 
its synthesis pp. However, very recent evidence from neutron diffraction indicates that 
below Tof the orbitals actually develop short-to-medium-range order in a trimerized 
state, although the associated strain field prevents the ordering from becoming long- 
range, and instead nanoscale domains of local orbital trimer order are formed |42j . 
Nevertheless, even local distortions favouring occupied d^z'^_r2 orbitals are intriguing, as 
we have seen above that the SE favours instead alternating symmetric and antisymmetric 
combinations of the basis orbitals, d^z^^r^ ± dx2_y2, for individual Ni^^-Ni^"*" pairs (so 
e.g. in a disordered state), and FO order of d^2_y2 orbitals, as obtained in the MF 
approximation. Actually, such an ordered phase with occupied dx2-y2 orbitals has been 
observed in large magnetic field |HI, which demonstrates that the state predicted by 
the SE model is energetically close to the actual ground state, and a quantum phase 
transition between different orbital states can be triggered by an applied magnetic field. 
The very fact that instead, in the absence of a magnetic field, local distortions of NiOg 
octahedra associated with occupied (ia;2_y2-type orbitals dominate at low temperature, 
we interpret as a strong indication that the orbital physics is in the end determined by 
the Jahn- Teller coupling to the lattice and not by the orbital SE. 

The ground state of LiNi02 remains the subject of intense debate. Several 
mechanisms for the spin liquid phase have been proposed so far: (i) quantum fluctuations 
on AF bonds ^3], (m) geometrical frustration of AF interactions in the S'[/(4) spin- 
orbital model on the triangular lattice [121 IHHI; iiH) decoupling of spin and orbital 
degrees of freedom because of the large difference, confirmed above, between orbital and 
spin-orbital SE coupling constants |2], {iv) disorder due to Ni^+ ions in the Li planes 
\^^. Interestingly, also in the SU{A) model both AF and FM interactions are allowed 
so at first sight the mechanism {ii) might appear similar to the present model of 
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section 0J However, the Cg orbital operators do not respect SU{2) symmetry, and thus 
the actual magnetic interactions which follow from the SU{A) model are quite different. 
For instance, in the S'f/(4) model one may assume perfect FO order, with (Tj ■ T,) = | 
for the orbital operators on each bond {ij), which would give AF spin interactions on 
all the bonds, while such a situation cannot be realized for eg orbitals. Also a state with 
the orbitals staggered for all bonds is not allowed on a triangular lattice. Therefore, we 
suggest that the role of geometrical frustration of magnetic interactions is overestimated 
in the S'f/(4) model and cannot explain the properties of the triangular Ni planes by 
itself. In contrast, all other mechanisms could contribute and stabilize the spin liquid 
state in a 3D lattice of LiNi02. We propose a new scenario below in which the above 
aspects (iii) and (iv) are supplemented by two important features of the orbital glass (TU] 
or trimer nanodomain state |l2] below the transition at Tof ~ 400 K: (v) the weakness 
of the (effective) SE interactions between Ni^"*" ions in adjacent Ni planes, and (vi) the 
Jahn- Teller coupling to the lattice, which induces local distortions favouring directional 
(i3^2_^2-like orbitals. 

We emphasize that the experimental studies on NaNi02 and LiNi02 have identified 
two competing effects which operate in triangular 2D Ni planes: on the one hand 
local Jahn- Teller distortions favouring occupied d^^'^-r'^ orbitals, produced by the on- 
site Jahn- Teller coupling, and on the other hand the SE, which would instead stabilize 
orbital order with occupied dx2_y2 orbitals, coexisting with FM interactions within the 
planes. Although Ni-Ni distances and Ni-O-Ni bond angles in the 0-Ni-O trilayers of 
the two compounds are similar, it may well be that the in-plane elastic couplings and 
orbital SE differ just sufficiently that the 2D orbital susceptibilities are qualitatively 
different. In particular, the competition between Jahn- Teller effect and SE may favour 
orbital disorder or complex multiple-sublattice orbital order when one interaction is 
only marginally stronger than the other one. This is in fact suggested by the behaviour 
of LiNi02 in magnetic field, mentioned above, and by the properties of the mixed 
compound Lio.3Nao.7Ni02 discussed below. Whether orbital and magnetic order 
stabilize in a 3D system, however, depends also crucially on the interplane interactions. 
Without any coupling between the planes, no long-range order could occur in the 2D 
planes anyway, in accordance with the Mermin- Wagner theorem. 

One thus expects that 3D orbital order can be induced only by a sufficiently strong 
effective interplane coupling, as apparently encountered in NaNi02, while the absence of 
orbital order observed in Lii_2;Nii+2^02 may be explained by a weaker effective interplane 
coupling in combination with an already weaker in-plane ordering tendency. Even 
without an explicit extension of the present microscopic spin-orbital model we can 
already suggest two mechanisms responsible for this weak interlayer coupling. First 
of all, the interplane SE is intrinsically weaker in LiNi02 than in NaNi02 because 
of the different covalency of the Ni-O-Li-O-Ni bonds as compared to the Ni-O-Na- 
0-Ni bonds, owing to the weaker hybridization between the oxygen 2p orbitals and 
the spatially less extended Li 2s orbitals than between the oxygen 2p orbitals and 
the bigger Na 3s orbitals. Although this difference is quantitative, it may lead to a 
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qualitatively different behaviour. A further point is then how the individual Ni^^-Ni^"*" 
orbital interactions between two neighbouring planes add up. For two orbital ordered 
planes, as in NaNi02, the interactions add up coherently, but if the orbitals preferentially 
avoid ordering in each Ni plane, e.g. because of a tendency towards domain formation, 
then in spite of the intrinsically identical spin-orbital interaction between Ni^+ ions in 
adjacent Ni planes, orbital (and magnetic) couplings add up incoherently. Second, it 
is well known that in the case of NaNi02 an almost ideal structure is obtained, while 
Lii_a.Nii+^02 is never stoichiometric and extra Ni ions occur in Li planes [0]. This creates 
frustration due to extrinsic randomness in the individual interplane interactions, and 
thus further weakens the effective interplane coupling. 

Finally, the absence of 3D magnetic order in LiNi02 is apparently similarly due 
to the weakness of the magnetic effective interplane coupling. While it is tempting to 
ascribe it instead to the presence of both AF and FM SE interactions in the triangular 
Ni planes, demonstrated above in the spin-orbital model, this is unlikely to be the 
leading effect stabilizing the spin liquid. In particular, in the orbital glass state the SE 
would be FM on by far the majority of the bonds in the Ni planes, as follows from 
equation [compare also figure |H] , with the precise fraction depending on the degree 
of disorder, while moreover the few AF interactions would be much weaker than the FM 
interactions. That the effective interplane magnetic coupling in LiNi02 is weak, is not 
only due to the interplane SE being intrinsically weak owing to the small covalency. In 
addition, in spite of the identical spin-orbital interaction between Ni^"*" ions in adjacent 
Ni planes, the absence of orbital long-range order in the Ni planes will make individual 
Ni3+_]\ji3+ spin-spin SE interactions add up incoherently. 

An additional extrinsic mechanism opposing magnetic order is provided by magnetic 
interactions due to the Ni^+ defects in the Li planes, which by themselves frustrate weak 
AF SE interactions by stronger FM interactions along Ni^^-Ni^^-Ni^"*" units ITT| . 
additionally enhanced due to the S = 1 spin states involved. As long as such defects 
do not occur, coherence in the magnetic interplane coupling is easier to obtain. Our 
scenario is not inconsistent with the recent observation of the decoupling of orbital 
and spin degrees of freedom in the mixed compound Lio.3Nao.7Ni02 which has a very 
different orbital state and yet a similar magnetic ground state as NaNi02 jH], with 
similar exchange constants j43j . 

8. Discussion and summary 

The above analysis demonstrates that the qualitative properties which follow from 
orbital and spin correlations within triangular Ni planes with 90° Ni-O-Ni bonds can 
be understood within a realistic spin-orbital SE model. This model demonstrates, in 
agreement with the SE model of Mostovoy and Khomskii [2 and with experiment |lUj . 
that the orbital SE is stronger by one order of magnitude than any other (pure spin 
or spin-orbital) interaction, because all magnetic dependence for the SE along Ni-O-Ni 
bonds originates from the singlet-triplet splitting of the oxygen 2p^ configuration, and 
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is therefore smaller by at least a factor Jp/Up ~ 0.1, where Up (Jp) is the interorbital 
Coulomb (Hund's exchange) interaction on oxygen. These parameters, as well as the 
effective hopping t and the charge transfer energy A, are of importance to establish 
quantitative consequences of the SE model, in particular on such magnetic properties 
as the effective exchange constants and the magnetic susceptibility. 

The orbital SE has interesting consequences for the orbital state and also for the 
magnetic interactions within the triangular Ni planes. First of all, for an individual 
(NiO)2 plaquette (figure Ej), a pair of identical orbitals, either both d3z2_r2 + dx2_y2 or 
both d3z^_r2 — d^2_y2, is favoured. Therefore, in MP approximation a symmetry broken 
state arises in the Ni planes, with FO order characterized by a single orbital angle 6 
[as defined in equation (j21)] along a particular direction. This is in contrast with the 
behaviour obtained for the cuprates or the manganites [HZj, where orbitals with 
angles 9 and —9 alternate on interlacing sub lattices. As a result, one expects domains 
with either \9) or \ — 9) orbitals, separated by twin-boundaries. Even though much 
weaker than the orbital SE, the spin-orbital SE plays a crucial role in the selection of 
the ground state. It breaks the U{1) symmetry generated in MF theory by the orbital 
SE and restores the original threefold symmetry with respect to the orbital angle 9. 
By this mechanism FO order with planar dx2_y2 orbitals in the triangular Ni planes is 
stabilized in MF theory. A characteristic feature of this state is the presence of both 
FM and AF spin interactions. 

After having a closer look at the magnetic properties of LiNi02, we suggest that 
this compound represents an interesting case of competition between the type of orbitals 
favoured locally by the Jahn- Teller effect and those favoured by the SE in an eg system 
with 90° Ni-O-Ni bonds, and that this is likely to induce orbital disorder. Remarkably, 
this case is reminiscent of the situation encountered in the t2g orbital model [201 1^ ■, and 
contradicts the experience and common knowledge from the Cg models in perovskites, 
where the SE along 180° bonds favours alternating orbitals and both mechanisms 
support each other US] • Whether this competition might lead to dynamical disorder 
certainly depends on the phonons, and we conclude that further study of elastic coupling 
and phonons is urgently needed. 

Summarizing, we presented and discussed the consequences of the 2D frustrated 
quantum spin-orbital model derived for the triangular Ni planes of two compounds with 
a very different behaviour: LiNi02 and NaNi02. This model emphasizes the importance 
of charge transfer (Goodenough) processes and presents a complete description of the 
superexchange interactions on the 90° Ni-O-Ni bonds in triangular Ni planes, so we 
are confident that it provides a solid starting point for future progress in the theory. 
The frustration in the orbital sector is not of Ising type, as suggested before jTU], but 
due to different orbitals being favoured for different bond directions. Providing a final 
answer concerning the origin of the different magnetic properties of LiNi02 and NaNi02 
requires an extension of the present microscopic model by a microscopic description of 
both the interlayer coupling and of the coupling between orbitals and the lattice. We 
identified two possible reasons why the interlayer (spin and orbital) coupling is so weak: 
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Figure Al. Example of superexchange processes (arrows in clockwise direction) of 
M-type along the XY direction. — After the first transition (top) either the triplet 
T or the orbital doublet D is reached at what is now a Ni^+ ion, while the second 
transition (bottom right) gives either the orbital doublet D or the singlet S at the 
other Ni2+ ion. The corresponding double excitation into the 2p'* configuration at 
oxygen is either to the triplet t or the singlet s. — Deexcitation occurs by the reverse 
transitions, either in the same or in reversed order (shown in the diagram). 

{{) an intrinsic effect due to the small covalency of interplane Ni-O-Li-O-Ni bonds, 
and (ii) an extrinsic effect due to Ni disorder in Lii_^Nii+^02 — both these effects 
are expected to play an important role and to be responsible for the absence of A-tjpe 
antiferromagnetic order in LiNi02. 
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Appendix A. Derivation of the spin-orbital Hamiltonian 

Here we give further details on the derivation of the spin-orbital Hamiltonian Ties, in 
particular regarding the case where Jh is finite. First we take a closer look at the four- 
step charge transfer sequences. As an example, let us consider (see figure lAljl a Ni-Ni 
pair in the XY direction in an M-type initial configuration with the electron on Ni-site 
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Figure A2. Diagrammatic representation of the superexchange processes shown in 
figure fSTl including inverted order of excitation and deexcitation processes. 



i (j) in the non-hopping (ij,2_^.2 orbital \y)i (hopping ds^^-r^ orbital \x)j) having spin up 
(down), i.e. an initial state |M, |J,), and the electrons being transferred from the oxygen 
having opposite spin. In an obvious notation the corresponding contribution to Hcs is 
[M, ti; itl VfVf Qjj. In the middle intermediate state the ion states T and D, t and 
s, and D and S occur with equal amplitude at the Ni^"*", 0° and Ni^"*" ions, respectively. 
Denoting the contribution from a doubly excited state by [Xt>y], with X,Y E {T, D, S} 
and V G {t, s}, the contribution corresponding to figure IXD is thus given by 

[M, ti; itl = U[TtD] + [TtS] + [DtD] + [DtS] + [TsD] + [TsS] + [DsD] + [DsS' 



where again a factor 4 has already been included. However, the four sequences that lead 
to a particular middle state upon inverting the order of the exciting and/or deexciting 
charge transfers, are now in general inequivalent, as illustrated by the tree diagram of 
figure IA21 because the energies of the first and third intermediate state may depend 
upon the sequence. The contributions [Xt>F] therefore need to be expanded as 



[XvY] 



(A.l) 



; {X\XvY\X) + {X\XvY\Y) + {Y\XvY\X) + {Y\XvY\Y) 

where we denote in {Zi\XvY\Z^) by Zi (Z^) the state reached after the first (third) 
charge transfer step. The explicit expression for this quantity is [compare equation 

m] 

{Z,\XvY\Z3) - 



Az,Az,i^x + AY + U, 



Az,Az,{^x + Ay) {Ax + Ay + U, 



Az,Az,iAx + Ay) 



(A.2) 
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where the second term in the first hne is the renormahzation correction discussed in 
section 121 Inserting ()A.2|) into equation ()A.1|) one now obtains 

«« Aa- + Ay U, 

SO that for example 

A U + J 

I'^'^l = 2(A2 - 4Jg)2 2A + C;, + J,' 

The next step consists in hsting all possible initial configurations and for each one 
the allowed spin arrangements of the two hopping electrons, and then working out, as 
above for [M, fj,; jf]; with the help of diagrams like those shown in figure IXTl which 
middle intermediate states contribute. The result is shown in table which gives the 
coefficients C {.]..) in 

[L, aa'; aa'l = ^ C(L, aa'; aa'|Xf F) [Xf F]. (A.6) 

X,Y,v 

The SE coupling constants {K'^'^} for spin triplet and spin singlet, appearing in equation 
(jnSl), are now obtained as follows. In each case (L = O, M or N) consider first the m = 1 
component of the triplet, i.e. the row(s) in table lAll in which there is an up-spin on 
both Ni3+ ions. Adding them gives the triplet contribution K^, while similarly adding 
the row(s) for which the two Ni^"*" spins are in up-down configuration gives Kj}, 

Kl = TT; aa'l = E E ^(L, TT; cra'lXvY) [XvY], (A.7) 

a,cr' a, a' X,Y,v 

Kl^ = Ti; aa'l = E E ^(L, Ti; aa'\XvY) [XvY]. (A.8) 

(7, a' a,a' X,Y,v 

As Kjj corresponds to the m = components and so is shared between singlet and 
triplet, one now obtains the singlet contribution from Kf^ = 2Kj} — K^. The result is 

= [DtD] + 2[DtS] + [StS], (A.9) 

Kl = [DsD] + 2[DsS] + [SsS], (A. 10) 

= 2[TtD] + [TsD] + 2[TtS] + [TsS] + [DtD] + [DtS], (A.ll) 

Ki^ = 3[TtD] + 3[TtS] + [DsD] + [DsS], (A. 12) 

= 7[TtT] + 2[TsT] + A[TtD] + 2[TsD] + [DtD], (A. 13) 

Kli = Q [TtT] + 3 [TsT] + 6 [TtD] + [DsD] . (A. 14) 

The above procedure (considering first the m = 1 component of the triplet) avoids the 
necessity of keeping track of the phases. Alternatively one may work out explicitly 
the states resulting after two perturbation steps from both |L, |i) and |L,i|), and 
from those, by taking their sum and difference, the states resulting from the m = 
components of the triplet and the singlet, and finally project the latter on all possible 
middle states. It is straightforward to verify that this gives the same results and 
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Table Al. CoefBcients C{L, aa' ; aa'lXvY) specifying the contributions to the 
superexchange Hamiltonian according to equation l)A.6(l . Here L, era' denotes the initial 
configuration and aa' the spins of the transferred electrons, with L e {0,M, N} and 
a, a', a, a' G {Tii}: ^^^d XvY denotes the middle intermediate state in the four-step 
charge transfer process, with X,Y E {T, Z?, 5*} and v E {t, s}. 
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equations ()A.9|) - ()A.14jl are reproduced. One should note that, while the parallel-spin 
initial state |L, ||) corresponds necessarily to the triplet and is therefore associated with 
FM coupling, it is not correct to associate the antiparallel-spin initial state |L, |J,) with 
AF coupling, since it projects on both the triplet and the singlet. 

The SE constants J^'^ occurring in equation (jH^ then follow from equations (jH^jl 
and (fXT9|l - |04|l . with the result 

j'o= + 2(D5)+ + (55) + , (A.15) 

J°= 3(rZ})+ + 3(TS)+ + (DD)++ {DS)+, (A.16) 

J° =9(TT)+ + 6(TZ})+ +{DDU, (A.17) 

Jg= iDD)_ + 2{DS)^ + iSS)_, (A.18) 

4= - (TD)_ - (T5)_ + + (DS)., (A.19) 

= (TT)^ - 2{TD)^ + (DD)^ , (A.20) 

where we have introduced the following abbreviations for the orbital exchange elements 
occurrmg m (|Xl5|l - (|XT7|l and the spin-orbital exchange elements occurring in ()A.18|) - 

(XF)+ = |[XtF] + i[XsF], (A.21) 
= [XtY] - [XsY]. (A.22) 

One readily verifies that equations ()A.15|) - ()A.20|) reduce to equations (PT |) - pi|) upon 
setting Jh = 0, whereupon all {XY)+ become identical as do all (XF)_. 
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From equations ()A.21|) - ()A.22|) or from the explicit expressions 

.yy. t' Ax + Ay {Up - iJpKAx + Ay) + Ul - Jl 



^ 2 AlAl- (Ax + Ay + U,y - 

it is obvious that the orbital exchange elements {(Xy),,.} are much larger than the 
spin-orbital exchange elements {{XY)^}, the ratio being 

|(Xy)_| 2J,(Ax + Ay) _ 4AJ, 

(f/p_ijp)(Ax + Ay) + ?72- J2 (2A + f/p)f/p' ^ • ^ 

where the second expression applies for Jp, J^i A,Up. It then follows from equations 
(|X?T5|) - (r07jl and (rO8|) - (rO0|) that the same holds for the orbital SE constant [or 
J!p, see equation ()49|l ] when compared with any of the spin or spin-orbital SE constants 
Jo-, or J^ [or Jr5, J^^., or J^g, see equation (jSHl)] — the same ratio ()A.25|) applies to 
Jts/ J'ti ill accordance with the analysis made above for Jh = [see equations (j43|l - (jl4j) ]. 

As for the relative size of the spin and spin-orbital constants with respect to one 
another, we observe that while Jq is the sum of four spin-orbital exchange elements 
[equation ()A.18j) ]. Jj^ is the sum of two differences of such elements, while is the 
difference of two differences of such elements, 

Ji = Ds + Dn, Ji = Dn- Dt, (A.26) 

with 

= (DX)_ - (rX)_ . (J- - J-)(J- + J- + ^)Cro.. (A.27) 
where Ctdx depends only weakly upon Jh and is given in good approximation by 

Obviously all Dx are positive since A^ < A^, and further Dt > Dd > Ds since 
At < Ad < As- It follows that J^ > and J^ < 0, and moreover, as clearly Ds + IDo 
is considerably larger than Dt-, one concludes that 

I^oI»4»I^nI- (A.29) 
Finally we arrive at the SE constants associated with the more physical interactions 
as occurring in TicS^o [equation P^jl ] and 7-^cfT,s [equation (|Hn|) ]. The orbital part is 
described by 

J. = i(J° - 2Jl + J°) = ^[9(TT)+ - 6(T5)+ + (55)+], (A.30) 

while the pure spin SE constant J^ and the two spin-orbital SE constants J^, and J^ are 
given by 

Ja = —-^("^o + 2Jm + -'n) 

= - ^[(TT)_ - 4(TD)_ - 2(T5)_ + 4{DD)_ + 4(D5)_ + (55)_], (A.31) 



Orbital and spin physics in LiNi02 and NaNi02 32 
= -i[(rT)_ + 2(TS)_ + (5S)_] 

= J, + 2Ds + 2Dd, (A.32) 
= - ^[-(TT)_ + 2(TD)_ + 2(D5)_ + (SS)-] 

= J„ + Ds + 2Dd - Dt. (A.33) 

From these expressions and the inequahty ()A.29|) above, the inequahty follows 
immediately. 
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